Draft version February 26, 2013 

Preprint typeset using I^'T^]X style eniulateapj v. 5/2/11 



CONFINED POPULATION III ENRICHMENT AND THE PROSPECTS FOR PROMPT 
SECOND-GENERATION STAR FORMATION 

Jeremy S. Ritter \ Chalence Safranek-Shrader \ Orly Gnat ^, Milos Milosavljevic \ and Volker Bromm ^ 

Draft version February 26, 2013 



(N 

o 

(N 



en 



O 
U 

6 

c3 



> 

in 

(N 

rn 
O 
(N 

> 

X 



ABSTRACT 

It is widely recognized that nucleosynthetic output of the first, Population HI supernovae was a 
catalyst defining the character of subsequent stellar generations. Most of the work on the earliest en- 
richment was carried out assuming that the first stars were extremely massive and that the associated 
supernovae were unusually energetic, enough to completely unbind the baryons in the host cosmic 
minihalo and disperse the synthesized metals into the intergalactic medium. Recent work, however, 
suggests that the first stars may in fact have been somewhat less massive, with a characteristic mass 
scale of a few tens of solar masses. We present a cosmological simulation following the transport of 
the metals synthesized in a Population HI supernova assuming that it had an energy of 10^""^ ergs, 
compatible with standard Type II supernovae. A young supernova remnant is inserted in the first 
star's relic H II region in the free expansion phase and is followed for 40Myr, all with the help of 
adaptive mesh refinement and Lagrangian tracer particle techniques. The supernova remnant remains 
partially trapped within the minihalo and the thin snowplow shell develops pronounced instability 
and fingering. Roughly half of the ejecta turn around and fall back toward the center of the halo, with 
1% of the ejecta reaching the center in ~ 30kyr and 10% in ^ lOMyr. The average metallicity of 
the combined returning ejecta and the pristine filaments feeding into the halo center from the cosmic 
web is '^ 0.001 — 0.01 Zq, but the two remain unmixed until accreting onto the central hydrostatic 
core that is unresolved at the end of the simulation. We conclude that if Population HI stars had less 
extreme masses, they promptly enriched the host minihalos with metals and triggered Population II 
star formation. 

Subject headings: cosmology: theory — galaxies: dwarf — galaxies: formation — hydrodynamics — 
ISM: structure — stars: formation — supernovae: general 



1. INTRODUCTION 

Early numerical investigations of the formation of the 
first stars in metal-free cosmic halos arrived at the con- 
clusion that the stars w ere isolated a nd rather massive, 
with masses > 100 A ^ (|Bromm et al. 1999, 2002; Abel 



stars ( |Bromm fc Loeb||2003[ [Smith et al]|2008[ ). 

Much of the theoretical work on the formation of the 
fir st stars and galaxies (e.g .. Wise fc Abel|2008a Wise et 



ar]|2012 Greif et aL]|2010 ), and on strategies tor detect 



et aL| 120001 |2002| |Gao et al.p()07| |0'5Eea ^TTTomla^ 



20071.^ This was a theoretically attractive possibility, 



mg metal-tree stellar populat ions in the observations of 

.,'Schaer er|2003l 
[Zackrisson et al. 



^laxies an d supernovae (e 



2009; Raiter et al. 2010 



Johns on et 
2UIi|, has 



given that such massive stars would have unique radia- 
tive, nucleosynthetic, and other signatures. For example, 
a metal-free 200 Mq star could emit 



postulated the massive character of metal-tree stars. Re- 
cent numerical investigations, however, have shown that 



ing photons per second (e.g., Bromm et al. 2001 Schaerer 
[2002 ), which would photoionize not only the star's par- 
ent cosmic minihalo, but a larger region of the primordial 
intergalactic medium (IGM). The star could explode as 
a ^ 10^^ ergs pair-instability supernova (PISN) , which 
could synthesize its own mass worth of metals. The 



angular-momentum-ind uced fragmentation (Stacy et al. 
2010|[Clark et al.ll2011|) and radiative termination of ac 



cretionj Hosokawa et a 



.|2011[[Stacyetal.|2012l see, also 



McKee fc ran||2008| during protosteiiar accretion could 



baryonic co ntent of the halo would be disrupted ( Bromm 



2003 ) and the metals would be dispersed far out- 
;hus enriching a r elatively large region of the sur- 



et al. 

side, 

rounding medium (e.g.. Wise fc A bel 2008b; Greif et al 
2010[[Maio et al .[[2010[ [Wise et al.||2012| ). While this 



enrichment would have only a minor efiect on the ulti- 
mate metal bu dget of the IGM ( Scannapieco et al. 20021 



Yoshida et al.[[2004| ), it would have the potential to raise 
the metallicity ot larger collapsing halos above a "floor" 
thought to be necessary for fragmentation into low-mass 
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limit the masses of the first stars to as low as a few tens 
of solar masses, and could yield initial protosteiiar mass 
functions co vering the entire range 0.1 — 10 M© (Greif 
et al. 20111. The final stellar mass range, at the end 
ot protosteiiar accretion, could be ~ 1 — 50 Mq with a 
characteristic mass of the order of 10 Mq. The cosmic 
imprint of such moderate mass first stars should be dif- 
ferent than that of the very massive ones: th e ionizing 
lumin osities could be < few x lO'^^ photons s~^ (Schaerer 
2002 1, supernova energies could be ~ 10^^ ergs, metal 
yields could be at most a few solar masses per super- 
nova, and the characteristic nucleosynthetic return could 
potentially be more co nsistent with currently available 
obser vational data (e.g., Tumlinson|2006 Vangioni et al. 



2011| ) 

The evolution of minihalos containing only moderate 
mass stars should differ in critical ways from that of ob- 



jects containing very massive Population III stars. With 
one-dime nsional hydrodynamic simulations, |Whalen et 
al. ( 2008 1 showed that the H II region of a moderate mass 
star at the center of the minihalo may remain confined 
within the minihalo. They calculated that the blastwave 
of a 10^^ ergs supernova in < 10^ Mq minihalos, does, 
after ^ 10 Myr, reach the virial radius of the halo, and 
concluded that this results in the "disruption" of the 
baryonic component of the halo. They also considered 
high-mass ~ lO'' Mq minihalos in which a strong pho- 
todis sociating background ha d prevented star formation 
(e.g., O'Shea fc Norman|2008 l, and found that there, the 
H II region may not form and the supernova blastwave 
may be confined at very small radii. They speculated 
that the subsequent prompt fallback of metal-enriched 
supernova ejecta into the center of the halo could lead 
to the formation of a globular-clustcr-likc stellar system, 
similar to the suggestion of Wada fc Venkatesan^ (2003 ) 
that a relatively small number ot supernovae exploding 
in a more massive, ~ 10® Mq cosmic halo could also lead 
to a fallback and prompt star cluster formation. This is 
a particularly interesting possibility, in view of the unex- 
plained origin of metal-poor globular clusters ( Brodie fc 
Strader||2006|: do they form in th e centers of their own 

e. 



low-mass dark m a tter halos 
blatt et al.|[T988l iBromm fc Clar 



2009 



Peebles 
"Ml 



ke 



.. ^^ [Boiey et~ 

Gritfen et al. 2010p, or do they forni ott-center in 



1984 Rosen- 



the merging and fr agmentation of more mas sive proto 



galactic disks (e.g., Kravtsov fc Gnedin||2005 and refer- 
ences therein)? 

The aim of this paper is to pursue implications of the 
hypothesis that the first stars had moderate masses — 
masses that are insufficient to produce extraordinarily 
extended H II regions and ultra-energetic supernovae. It 
will be necessary to understand the entire evolution of 
supernova ejecta, from the point of injection in the su- 
pernova explosion until some fraction of the ejecta has 
become incorporated in new stars. Existing treatments 
of the initial phases, including metal injection in a su- 
pernova, the subsequent excursion into the environment, 
and a possible fallback in t he host halo, have entailed 
one-dimensional simu lations ( Kitayama fc Yoshida|200"5" 



Whalen et al. 20081 or relatively low-resolution three- 
dimensional, grid- based or smoothed particle hy dro dy- 
namic simulations tNakasato & Shigeyama 2000 ; Bromm 



et al.|2003[[Greif et al.|2007[|Wise fc Abel|2008b|) . These 
treatments point to either a prompt and widespread 
dispersal of the supernova ejecta into the primordial 
IGM — assuming that the supernovae were extraordinar- 
ily energetic — or, in the case of particularly high-mass 
halos hosting moderate mass stars, to a possible suppres- 
sion of H II region breakout and the subsequent trapping 
and prompt fallback of the ejecta. 

It is important to realize, however, that the fate of 
the supernova ejecta may be influenced by dynamical 
processes operating on relatively small spatial scales and 
strongly violating spherical symmetry. In addition to the 
likely intrinsic asphericity of the explosion itself, one such 
process is the Rayleigh- Taylor instability of the decelerat- 
ing contact discontinuity after the ejecta have started to 
interact with the ambient medium. The ambient medium 
into which the supernova expands is bound to be inho- 
mogeneous, especially if the luminosity of the progenitor 
stars and other hot stars that have formed with it is in- 



sufficient to completely photoionize the densest clumps 
remaining from the progenitor gas cloud. On the some- 
what larger ~ 100 pc scales of the host minihalo, a pro- 
nounced departure from spherical symmetry will arise 
from the filamentary geometry of the cosmic web. Such 
asymmetry can interact with the expanding blastwave. 
Further, as an aging supernova remnant forms a thin, 
momentum-conserving snowplow shell, the decelerating 
shell is itself unstable to the Rayleigh- Taylor instability 
and forms fingers. Finally, while the blastwave insta- 
bilities lead to partial mutual entrainment of the ejecta 
and the shocked primordial gas, any ejecta-enriched gas 
falling back toward the center of the halo, where it may 
engage in a new round of star formation, may not nec- 
essarily be chemically homogeneous (full chemical mix- 
ing requires at least a few turbulent eddy turnover times 
which may not be available in the blastwave). 

In view of these complexities, here we report the first 
results from a program to study the detailed transport 
of the nucleosynthetic output of the first stars, from the 
point of injection in supernovae, until at least a fraction 
of the new metals have recondensed in a descendent cos- 
mic halo where they may contribute to the formation of 
second-generation, metal-enriched stars. The character 
of the new star formation in the descendent halo must 
depend on the rate of gas inflow Minf into the nucleus 
of the halo, where the gas can rapidly cool and fragment 
into stars. The cooling rate, and the compressibility of 
the turbulence that may be present, are sensitive to the 
metallicity of the inflowing gas Zjnf and to the degree of 
metallicity homogeneity in the gas. We present the re- 
sults of a cosmological simulation specifically designed to 
measure these parameters and thus shed light on the con- 
ditions for the formation of metal-enriched stars in the 
aftermath of a Population III supernova, assuming that 
it exploded with only 10^^ ergs of ene rgy. The simulation 



shares some similarities with that of Boley et al. ( 2009 1 
and its setup has many common elements with tnat of 
Wise et al. (2012). The latter work followed the evolu- 
tion ot a cosmological region of the same size as ours but 
to much lower redshifts, with a focus on the statistics 
of star formation and metal enrichment across a large 
number of star-forming halos. Our focus, instead, is on 
the detailed, smaller-scale character of metal transport 
in the aftermath of a single supernova. 

The paper is organized as follows. In Section [2] we de- 
scribe our numerical algorithm. In Section [3J we present 
our results. In Section HI we discuss implications for 
the formation of metal-enriched star clusters in the early 
universe, and in Section [Sl we present our main conclu- 
sions. We adopt the cosmological parameters consistent 
with the Wilkinson Microwave Anisotropy Probe seven- 
year data (Komatsu et al. 2011). We explicitly indicate 
whether comovmg or physical length units are utilized; 
all other quantities, such as density, are expressed in 
physical units. 

2. SIMULATING A SUPERNOVA IN A MINIHALO 

To investigate the fate of the first metals injected into 
the universe by a moderate mass Population III star, we 
carry out a cosmological simulation large enough not only 
to permit the formation of such a star, but also to allow 
the subsequent growth of the star's host halo, by accre- 
tion and merging, into a larger object that is the poten- 



tial site of second-generation star formation. The metals 
must be inserted early on, before the freely expanding su- 
pernova remnant has started to interact with the primor- 
dial circumstellar medium; this is possible thanks to our 
ability to drastically increase the local mesh resolution at 
a specific time in the simulation — ^just prior to supernova 
explosion — and then gradually degrade the resolution as 
the remnant expands. To reproduce the remnant's dy- 
namics correctly, during the progenitor star's life on the 
main sequence, we calculate the transfer of the star's 
ionizing radiation and the response of the circumstel- 
lar medium by the associated photoheating in the H II 
region. The simulation, which we proceed to describe, 
is designed to resolve the rudimentary time and length 
scales required to assess the impact of the supernova 
on the host halo and the nature of the metal-enriched 
flows in the aftermath of the supernova: how quickly do 
baryons return into the center of the halo (or one of its 
descendent halos), and how do the baryonic accretion 
rate and the metallicity of the accreted baryons evolve? 
The descrip tion of the simulation is organized as fol- 
lows. Section |2.1| describes the procedure by which the 
initial conditions were generated and the gravitational 
potential computation. Section |2.2| describes the chem- 
ical and thermodynamic cal cula tions, with and without 
an ionizing source. Section 2.3 describes the strategies 
we employed to control the mesh refinement in the sim- 
ulations. Section [2 .4| describes our algorithm for transfer 
of i onizing radiation from the first star. Finally, Section 
[23] describes the insertion of the supernova remnant and 
the subsequent tracking of the metal transport in the 
simulation. 

2.1. Cosmological Initial Conditions and Gravity 
We initialize our simulation in a periodic cosmological 



grid. 



box of volume 1 Mpc (comoving) at redshift 



146. 



The initial conditions are generated with the multi- 
scale cosm ological initial conditions package GRAPHIC2 
( Bertschinger 2001) with two levels of nested refinement 
on top of the 128'^ base grid. The highest refined region 
had an effective resolution of 512'^, corresponding to a 
dark matter particle mass of 230 M© . While we will 
focus on the dynamics of an initially ~ 10^ Mq minihalo 
at redshift z ~ 20, the region of high refinement was 
centered on the overdensity that first collapses to form 
a more massive ~ 10^ Mq halo at a lower redshift. The 
time integration was carried out with the multipurpose 
astroph ysical adaptive mesh refinement (AMR) code 



FLASH (Fryxell et al. 2000), version 3.3, em ploying 
the direct multigrid Foisson gravity solver of [Ricker 
( [2008 ) . Th e AM R refinement criteria, which we describe 
in Section |2.3| below, are set to ensure that relevant 
baryonic processes are sufficiently well resolved in the 
simulation. In collapsing regions where baryons dom- 
inate the gravitational potential, the AMR resolution 
length quickly becomes smaller than the dark matter 
particle separation. There, we employ a mass- and 
momentum-conserving quadratic-kernel-based so f tening 
procedure developed in Safranek-Shrader et al. (20121 
which is applied to the dark matter density variable 
during mapping onto the computational mesh, to render 
the dark matter particle contribution to the gravita- 
tional potential smooth on the scale of the local AMR 



2.2. Chemistry and Cooling 



Frior to the initialization of the H II region, we inte- 
grate the full nonequilibrium chemical network for hydro- 
gen, helium, deuterium, and their chemical derivatives. 
The chemical and radiative cooling updates, which in- 
clude, among other processes, the cooling or heating by 
the cosmic microwave background (CMB), are operator- 
split from the hydrodynamic update and are subcycled 
within each cell. The chemical reaction rates are the 
same as we have summarized in ISafranek-Shrader et aL] 
(20101, except that the chemical and radiative coohng 
update is isochoric in the present simulation. The net- 
work is integrated with the Bulirsch-St oer-type semi- 
implic it extrapolation method of Bader fc Deuflhard| 
([1983). W e test the therniodyna mic evolution of the gas 
in , Safranek-Shrader et al. (2012 1. 

After turning o n th e ionizing source, which we de- 
scribe in Section 2.4 below, we continue to integrate 
the nonequilibrium chemical network in the neutral gas, 
but in the ionized gas, we adopt a scheme in which the 
gas temperature and chemical abundances relax expo- 
nentially toward a target photoionization equilibrium. 
The equilibrium gas temperature and chemical abun- 
dances are functions of the photoionization parameter 
^ = AnF/nn only, where F is the radiative energy flux 
and nil is the number density of hydrogen nuclei and are 
precomputed with the same chemical network that we 
integrate in the neutral gas. The relaxation time scale is 
chosen to be the true, nonequilibrium thermal time scale 
ithcrm = {dhiT/dt)~^ , defined by the instantaneous local 
radiative cooling and photoionization heating rates. 

After inserting the supernova remnant, the gas can 
contain an arbitrary admixture of metals in addition to 
the primordial elements. The radiative point source is 
removed, thus discontinuing the photoionization equi- 
librium calculation. We revert to an integration of 
the nonequilibrium chemical network for the primordial 
species and calculate the ionization state of the met- 
als and their contribution to the cooling by interpola- 
tion from precomputed tables. At temperatures above 
8000 K, the collisional ionization equilibrium state of 
each of the metal species, including C, N, O, Ne, Mg, 
Si, S, and Fe, is calculated separate ly us ing methods de- 
scribed inlGnat & Sternberg[ ( 2007| and [Gnat & Ferland 



(20111. Below HUOUK but aBwe 200K, we ignore the 



metallic contribution to the free electron abundance and 
tabulate a single density-independent cooling function 
(such that the volumetric radiative loss rate scales with 
the square of the density) , appropriate to a gas at rela- 
tively low densities 5 , lO'^cm^^, using the code CLOUDY 



(Ferland et al. 



minimum met a 



1998[ ). The tabulated data extend to a 
licity of 10"'^ Zq; below this metallicity. 



we set the metal cooling rate in this range of tempera- 
tures to zero. Below 200 K, we disable cooling by metal 
lines altogether; however, because low-temperature ther- 
modynamics is not the focus of this work, this has no 
effect on the forthcoming results. The treatment of pri- 
mordial species takes into account the electrons provided 
by the metals and vice versa, but we do not include colli- 
sional charge exchange reactions between the primordial 



species and the metal species. We assume that the metals 
have solar abundance ratios. This is an arbitrary, unnat- 
ural choice given that our metals are being produced by a 
core-collapse supernova, but this choice should have only 
a very minor effect on the thermodynamic evolution of 
the supernova remnant. 

2.3. Adaptive Mesh Refinement Strategies 

The simulation is set to adjust the FLASH refinement 
level £ based on the local baryon density pb to ensure 
that denser regions are more highly refined. The grid 
separation (cell size) is related to the refinement level via 
Ax = 2~^~^^~^L, where L is the size of the computational 
boxjj As gravitational collapse drives the local fluid den- 
sity pb to increasingly higher densities, the local refine- 
ment level is adjusted to satisfy pb < 3pb 2^^^"'"''^^^^~^''"'''' 
at all times, where ^base = 5 is the refinement level corre- 
sponding to the base grid of size 128 = 2^'"=''"'"^, and con- 



sistent with Safranek-Shrader et al. (2012 1, we choose the 
scaling parameter to equal = — U.3. We insert the star 
particle when the hydrogen number density first exceeds 
1000 cm^^; this happens at redshift z « 19.7 and refine- 
ment level £ = 14, corresponding to cell size Ax sa 0.74 pc 
(physical) . 

At the end of the star's life and just prior to inserting 
the young supernova remnant, we force the simulation to 
refine further to level ^ = 21 at the location of the star 
particle, corresponding to cell size Aa; sa 0.006 pc (physi- 
cal) . This allows us to insert the supernova remnant well 
within the free-expansion phase, before the ejecta have 
started to interact with the circumstellar environment. 
Then, we disable automatic derefinement, and allow the 
simulation to refine additional computational cells, up 
to the maximum level ^max initially equal to 21, based 
on the standard second derivative test in FLASH, but in 
this case, the second derivative test is computed using 
the cell metallicity. This ensures that the entire region 
enriched with metals will be resolved at the refinement 
level ^max- As the supernova remnant expands, we grad- 
ually and cautiously lower ^max, and thus coarsen the 
resolution of the metal-enriched region, to allow us to 
make ever longer time steps and integrate the simulation 
for 40 Myr past the insertion of the supernova. 

2.4. The First Star and Ionizing Radiation Transfer 

When the hydrogen number density first exceeds 
1000 cm^"^ near the center of the collapsing minihalo, 
we insert a collisionless star particlen The star is free 
to move in the combined gravitational potential of the 
dark matter and baryons but does not contribute to 
the potential, i.e., it is treated as massless. We assume 
that it is an isotropic source of hydrogen-ionizing radi- 
ation and treat the radiation as monochromatic with 
photon energy hv = 16 eV and photon emission rate 
(5(H) = 6 X lO"*^ photons s""'^. This ionizing photon rate 
is about three times as high as has been estimated for a 
40 Mq metal-free star without mass loss, Q{li,4:0 Mq) « 

^ The factor of 2""^ arises from the subdivision in FLASH of each 
computational block relative to which the refinement level is de- 
fined into 8^ cells. 

^ This threshold density for star formation is ad hoc, but con- 
sistent with the choice in, e.g.,|Wise et al.|(|2012||. 
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Figure 1. The halo mass M200 and radius -R200 as a function 
of redshift, over the course of 50 Myr. The dotted line marks the 
formation of the star, and dashed line the supernova explosion. 
The steps are an artifact of the method we use to estimate iJ200. 



2 X 10**^ photons s^i ( |Schaerer||2002| . The ionizing pho- 
ton energy is calibrated to the latter lower luminosity 
star. We deliberately overestimate the ionizing luminos- 
ity to allow for the presence of additional, lower-mass, 
longer-lived stars also producing additional ionizing pho- 
tons in a compact cluster around the primary 40 Mq star. 
Also, the artificially high ionizing luminosity allows the 
H II region to break out more easily at the moderate spa- 
tial resolution at which we insert the star particle and 
perform ionizing radiation transfer. 

To carry out ionizing radiation transfer, we center a 
spherical coordinate system on the current position of the 
star particle, and pixel ize the angular coo rdinate using 
the HEALPix scheme of |G6rski et al.| ( |2005| ) with A^^pix = 
3072 pixels. Within each angular pixel, we split the ra- 
dial coordinate in the range 4 x 10^° cm < r < lO^'' cm 
(comoving), corresponding to 0.7pc < r < 17kpc (phys- 
ical), into iVrad = 103 logarithmic radial bins. For each 
volume element Epj, defined by an angular HEALPix pixel 
p and a radial bin 6, and for each mesh cell Ec defined 
by the cell index c, we compute the volume of the inter- 
section of the two elements V^,&;c = vol(i?p.b n Ec) > 0. 
In general, only a fraction of the cell lies within Epj,. In 
this case, we recursively split the cell into sub-cells via an 
octree subdivision procedure and compute the contribu- 
tion of the subcells lying within Ep,i, to the intersection 
volume; the recursion is continued until a desired accu- 
racy is achieved. Armed with the intersection volumes, 
we compute the average case-S hydrogen recombination 
rate n^cc,p,b evaluated under the assumption of full ion- 
ization (as in the standard Stromgren calculation) within 
each Ep^b via 



\cc,p,b 






arcc,B(T'c)"'c^P,&;c 



(1) 



where Vp^b = vo\{Ep^b) ~ J2c ^pMc: Tc and n^ are, re- 
spectively the temperature and density of hydrogen nu- 
clei within the cell, and arec,B(T) is the hydrogen re- 
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Figure 2. Density-temperature phase plot for the gas within 1 kpc (physical) from the center of the minihalo at redshifts, from left to 
right, z = 19.7 (right before the insertion of the star), z = 19.5 (in the Sedov- Taylor phase of the supernova, 8kyr after the explosion), 
and at z = 17 (end of the simulation, 40Myr after the explosion). By the end of the simulation, the gas heated by photoionization and 
the supernova has cooled and the halo has resumed central gravitational collapse, albeit at four times higher halo mass and now a nonzero 
metallicity. The dashed line shows the CMB temperature. 



combination coefficient. Within each bin, we compute 
the approximate integral of 47rr^nrcc,p('") along the ra- 
dial direction (now replacing the bin index b with the 
continuous radial coordinate r) and determine the pixel- 
specific Stromgren radius Rs,p that solves the equation 



Q(H) = 



Rs,. 



Aiir h 



rcc,p 



{r)dr. 



(2) 



We assume that no ionizing radiation penetrates to radii 
r > Rs,p- At smaller radii, we compute the local ionizing 
photon number flux (in the units of photons cm~^ s^^) 
using 



frir) 



1 



Q(H) 



'i'n-r^hrcc,p{r')dr' 



(3) 



The local ionizing flux is then used to compute the local 
photoionization parameter for cell c lying within the pixel 
p with a center at a radius r^ via S.p{r) = 4:T:hi>fp{rc)/nc. 
Given a value of the photoionization parameter, the pho- 
toionization equilibrium temperature of the gas and the 
chemical abundances within the cell are evaluated via 
table lookup (see Section 2.2). The computation is fully 



parallel and is repeated at every hydrodynamic time step. 
We do not compute the radiation pressure from Thomson 
scattering and photoionizations. This radiative transfer 
scheme is not explicitly photon conserving and is not de- 
signed to correctly reproduce the ionization front kine- 
matics, but it is sufficiently accurate in the present con- 
text where the sole interest is in the hydrodynamic re- 
sponse of the gas to photoionization heating and its im- 
pact on supernova blastwave dynamics. 

2.5. Supernova and Metal Transport 

We insert the supernova in the free-expansion phase, 
when the remnant is about t — igN = 35 yr old. The 
ejecta are cold and have a total mass AfsN — 40 Mq and 
a total kinetic energy of i^sN — 10^^ ergs. The metal 
yield is assumed to be ZgN — 0.1 and uniformly pre- 
mixed in the ejecta, implying 4 Mq of metals. The ve- 
locity of the ejecta is linear in the distance from the cen- 
ter of the explosion and the density is uniform. The 
initial radius of the remnant is chosen such that its ra- 
dius is 10% of the radius of the remnant at which the 
swept up mass equals the ejecta mass, or in this case, 
-RsN = 0.075 pc. The simulation carries out advection 



of the absolute metallicity Z by treating it as a pas- 
sive mass scalar quantity. This advection is subject to 
undesirable numerical diffusion, which is unavoidable in 
Euler ian solvers of this type (see, e.g., iPlewa & Miiller 



19991. Therefore, in parallel with passive scalar advec 



tion, we also insert a number iVpart = 5 x 10 of pas- 
sive Lagrangian tracer particles, originally uniformly dis- 
tributed in the ejecta. The simulation advances particle 
positions by integrating the first-order ordinary differ- 
ential equation dH-i/dt = v(Xi), where X; is the loca- 
tion of particle i, and v(Xi) is the baryon fluid velocity, 
quadratically interpolated from the computational grid, 
at the location of the particle. We do not take into ac- 
count the possibility that the supernova may be a source 
of molecules and dust. Both are undoubtedly produced 
in the free-expansion phase, but their survival of the re- 
verse shock hinges on small-scale dumpiness of the ejecta 
and on ot her complex, incompletely-unde rstood physics 



(see, e.g., Cherchneff &: Dwek 



2009 



2010) 



3. RESULTS 
3.1. Thermodynamical Properties 

In Figure[l} we show the redshift evolution of the radius 
i?2oo within which the average halo density equals 200 
times the critical density, as well as the corresponding 
enclosed halo mass M2oo- The star particle is inserted 
at z « 19.7 and the supernova at z « 19.5. The halo 
radius doubles from i?200 ~ 160 pc to 300 pc in the course 
of the simulation, and the halo mass quadruples from 
M200 « 10^ Mq at star insertion to 4 x 10'' Mq at the 
end of the simulation. 

Figure [2] shows the density-temperature phase diagram 
for gas within 1 kpc (physical) from the center of the 
halo, immediately prior to the formation of the star par- 
ticle, after the supernova explosion, and finally at the 
end of the simulation when the halo mass has quadru- 
pled. Comparison of the thermodynamic profile before 
stellar insertion and at the end of the simulation illus- 
trates how the halo relatively quickly "recovers" from the 
radiative and kinetic impact of the first star and resumes 
central gravitational collapse. Prior to the photoioniza- 
tion by the star, the archetypical thermodynamic struc- 
ture of a minihalo in which gas cools by rovibrational 
transitions of molecular hydrogen is evident. Photoion- 
ization and the supernova blastwave temporarily disrupt 
this structure. The photoionization raises the tempera- 



1000.0 



o 

Cl 




0"" 10' 
t (Myr) 

Figure 3. Characteristic radii corresponding to Lagrangian mass coordinates enclosing baryonic masses of, from top to bottom, 
(2.1, 1.0, 0.5, • • •) X 10® AIq (left panel) and ejecta mass fractions, as tracked by passive tracer particles, of (100%, 50%, 25%, • • ■) (right 
panel), both as a function of time since the explosion. The dotted line shows the virial radius R200 of the host halo. The baryonic radii 
are computed from the gravitational potential minimum of the halo. Supernova-driven expansion is evident in the mass coordinates inside 
the virial radius, but all mass coordinates return to infall owing in part to a persistent filamentary streaming of fresh, cold baryons from 
the cosmic web into the halo center. For ejecta radii, at times t < 0.02 Myr, the location of the progenitor star is taken as the coordinate 
origin; at later times, the origin is the location of the gravitational potential minimum; the change of origin is visible as a discontinuity in 
the curves. By the end of the simulation, a half of the ejecta has turned around and is falling back toward the center of the halo. 



ture of the photoionized gas to « 10'* K, while the explo- 
sion introduces hot gas that cools both adiabatically and 
radiatively. Because of the low energy and compactness 
of the remnant, the cooling by inverse Compton scat- 
tering of the CMB is relatively unimportant, different 
from the case of ultra-energetic Population III super- 
novae, wh ere this cooling cha nnel dominates the early 
evolution ( Bromm et al.|[2003 ). 

It is worth noting that an examination of the three- 
dimensional structure of the thermodynamic state of the 
gas shows that the dense, cold gas inside the filamentary 
structure feeding from the nearby cosmic web into the 
center of the minihalo is self-shielding and remains un- 
affected by photoionizationj^ The blastwave also largely 
avoids these dense filaments, and they remain a reposi- 
tory of low-entropy gas. This is easily seen in Figure HI 
which shows mass-weighted projected density, tempera- 
ture, and metallicity 8.5 Myr after the explosion. Outside 
of the filaments, after 40 Myr, the effects of photoioniza- 
tion and the blastwave have largely been erased; the ther- 
mal phase structure within the halo is similar to that of 
a higher-mass unperturbed minihalo. The appearance of 
warmer (< 10^ K) gas at low densities refiects the higher 
virial mass of the halo at the end of the simulation and 
the residual entropy left by the fossil H II region. The 
presence of metals at local metallicities Z > 10~^ and an 
enhanced abundance of hydrogen-deuteride now facili- 
tate more rapid cooling toward low temperatures. The 
enhanced cooling is now especially apparent at high den- 
sities riH ^ 100 cm""^ that are found in the central region. 
The dense gas is insufficiently well resolved at the end of 
the simulation; the resolution limit hinders compression 

^ The simulation does not properly resolve the length scales in- 
volved in the photoevaporation of dense cloud cores, thus we are 
not sure that the simulation has converged in this respect. 



of the fluid to densities tih ^ 1000 cm^"^. Nevertheless, 
it is clear that the dense gas resides at temperatures a 
factor of a few lower than the gas at the same densities 
prior to the insertion of the ionizing source. We do not 
observe the anticipated gas cooling to the CMB temper- 
ature, either through metal fine cooling or the coolin g 
by hydrogen-deuteride (e.g., [Johnson fc Bromm| 2006), 
because we enforced an artificial temperature floor at 
200 K. 

3.2. Blastwave Evolution and Fallback 

An insight in the rate of recovery of the halo can be ob- 
tained from Figure [3] which shows the Lagrangian radii 
that track the radial location of a grid of mass coordi- 
nates of the spherically averaged baryon and supernova 
ejecta mass as a function of time from the explosion. 
It is worth keeping in mind that the underlying flow is 
complex and contains simultaneous inflows and outflows 
that are averaged over; each of the radii reflects only the 
net displacement of the mass coordinate. The baryonic 
Lagrangian radii show expansion, driven by the super- 
nova blastwave, starting at w 20kyr at the innermost ra- 
dius shown, and appearing at progressively larger mass 
coordinates until outward-directed motion becomes evi- 
dent near the virial radius at « 2 Myr after the explo- 
sion. Each of the radii corresponding to mass coordi- 
nates inside the halo, however, shows a reversal from a 
net outflow to a net inflow. The reversal commences 
at « 400 kyr at the innermost radius shown, and reaches 
progressively larger radii until ~ 20 Myr when the flow at 
the virial radius turns around into an inflow. Outside of 
the virial radius, cosmic infall persists largely unabated 
throughout the remnant's history. 

An examination of the evolution of the actual three- 
dimensional structure of the baryonic flow (see also the 




Figure 4. Baryon mass-weighted projection of hydrogen density, temperature, and fluid-based metallicity in units of the solar metallicity 
8.5 Myr after the explosion. The cirele shows the virial radius of the halo, which is here 180 pc (physical). 



projection in Figure HI) reveals that the prompt reversal 
can only in part be attributed to a fallback of the mass 
that the blastwave has swept up and to which it has 
transferred momentum. Instead, the filamentary inflow 
of dense, cold, molecular gas from the nearby cosmic web 
directly into the center of the minihalo contributes a large 
fraction of the mass inflow rate. The baryonic filaments 
are seen diagonally from top left to bottom right in Fig- 
ure [Sl which shows a slice of the density field. Because 
the filaments survive photoionization and the explosion 
largely intact, they can sustain a net inflow into the cen- 
ter of the halo even when a majority of the ejecta and 
the swept up mass are moving outward. 

The Lagrangian radii for the supernova ejecta in Fig- 
ure [31 right panel, show first signs of interaction of the 
ejecta with the circumstellar medium around 1 — lOkyr 
after the explosion. What follows is a similar picture 
of expansion eventually reversing into fallback. Frac- 
tions of (3.125%, 6.25%, 12.5%, 25%, 50%) of the ejecta 
mass turn around and start falling back at around 
(0.2, 1.0, 4.0, 8.0, 3.0) Myr after the explosion. The cor- 
responding turnaround radii are (3, 8, 30, 60, 200) pc. 
The outer somewhat less than ^ 50% of the ejecta passes 
the virial radius and continues to travel outward until 
the end of the simulation. The net outflow of ejecta at 
T ^ ^200 can be contrasted with the net inflow of baryons 
at the same radii, underscoring the three-dimensional na- 
ture of the flow and an inadequacy of one-dimensional 
integrations in treating the long-term dynamics of su- 
pernova remnants in cosmic minihalos. The spherical 
symmetry implicit in one-dimensional integrations forces 
the ejecta and the swept up circumstellar medium into a 
spherical thin shell at late times. It does not allow for the 
presence of the dense, low-entropy clouds that enter the 
halo by infall from the fllaments of the cosmic web. The 
expanding shell may pass and leave behind such clouds. 
The spherical symmetry also does not allow radial seg- 
regation of the ejecta as a result of the instability of the 
decelerating thin shell. 

Rayleigh- Taylor flngering is evident in Figure |4] at 
8.5 Myr after the explosion. The long term evolution 



of the fragmented shell is seen in Figure [5) left panel, 
which shows the ejecta tracer particles in projection, su- 
perimposed on a slice of the density field through the 
center of the halo. It is also clear in these figures that 
the blastwave has expanded biconically perpendicular to 
the sheet-like baryonic overdensity deriving from the cos- 
mic web; this pristine gas, unpolluted by the ejecta, is 
still able to stream into the halo center. In the right 
panel of Figure [5j showing the central 360 pc, pristine 
baryonic streams arrive diagonally into the central few 
tens of parsecs, where they form a two-armed spiral and 
join a self-gravitating central core. At this point in the 
simulation, the spatial resolution is insufficient to resolve 
the internal structure of the core. 

Strikingly, we observe filamentary accretion of ejecta- 
enriched gas perpendicular to the pristine sheet. This 
is the return of the ejecta from the supernova remnant 
back into the center of the halo. A closer examination of 
the time evolution of the ejecta kinematics shows that 
in the snowplow phase, the ejecta are compressed into 
a thin shell, as expected. The thickness of the shell is 
itself not well resolved, thus we are not in the position to 
explore its s usceptibility to direct gravitational fragmen- 



tation (e.K.,|Salvaterra et al. i2004; Machida et al.||2005, 
Chiaki et al.||2012[ ). TIie"Raylcigh-'i'aylor instabflity 



buckles and redistributes momentum in the shell. The 
highest momentum fragments continue semi-ballistically 
outwa. r d, co nsistent with the observation by |Whalen| 



et al. (2008J) that the blastwave of even a moderate 
j^QOi gj.gg explosion has enough momentum to escape 
a standard minihalo. The low-momentum fragments, 
however, do not have sufficient momentum to escape; 
they turn around, fall back, and eventually join the 
central unresolved core. The falling ejecta make thin 
multiply-folded streams and filaments, indicating what 
seems to be a certain lack of mixing with the halo gas. 
The pristine gas accreting through the filaments of 
the cosmic web, and the ejecta-enriched gas accreting 
from the supernova remnant, meet in the central core. 
If the core is turbulent — poor resolution at this stage 
of the simulation would not allow us to detect such 
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Figure 5. Metal dispersal and fallback at the end of the simulation, 40Myr after the supernova explosion, at z = 17. The left panel 
is 1.1 kpc (physical) wide and centered on the gravitational potential minimum; the right panel is a 360 pc detail. The black points are 
the metal particles in projection, while the color shows a slice through the physical hydrogen density njj- Metal-bearing Rayleigh- Taylor 
fingers have a positive radial velocity and have breached the virial radius of the halo, while a majority of the metal mass is falling back 
into the halo center and remains incompletely mixed with the primordial gas. Density in the central unresolved core exceeds coverage of 
the color scale and is ~ 10* cm~^. 



turbulence — turbulent mixing of the ejecta and some 
swept up pristine gas with a much larger mass of pristine 
gas would ensue. 

3.3. Metal Enrichment 

In Figure [6] we show the net rates of baryon and metal 
outflow from or inflow into the central 20 pc of the halo. 
Some baryonic inflow is already evident early on, after 
50kyr from the explosion, and after « 1.5 Myr, baryons 
transition to permanent inflow. After w 5 Myr, the in- 
flow rate settles at Mb « 0.002 Moyr"^ The metal 
flow transitions from outflow to inflow at ~ 4 Myr af- 
ter the explosion, and fluctuates in the range Mz ~ 
(0.5 — 5) X 10~^ Mq yr~^. The average asymptotic abso- 
lute metallicity of the inflowing material is 



^inflow 



Mz 
Mb 



10" 



"oT 



0.005 Z, 







"oT 



.(4) 
where we are referring to the total metal mass ignoring 

fractional abundances of individual metal species. 

We proceed to examine the evolution of the metal- 
mass-weighted metallicity probability density function 
(PDF). The metaUicity PDF is computed in three differ- 
ent ways. In the first, "fluid-fluid" method, the metallic- 
ity is obtained by summing the amplitudes of the passive 
scalars defining metal abundances on the computational 
grid. The metal mass at metallicity Z is then calcu- 
lated by computing a volume integral of Zp over the cells 
with metallicities in a narrow logarithmic bin containing 
Z. The shortcoming of this method is that it yields, as 
a consequence of a well-known shortcoming of passive 
scalar transport schemes, nonvanishing, small metallici- 
ties even in many cells that ejecta could not have reached 
on hydrodynamic groundsr] With this in mind, the first 

* The spurious diffusion can easily be understood by considering 



method is more accurate where the metallicity is near its 
peak value, and is highly inaccurate elsewhere. Our re- 
maining two methods for computing the metallicity PDF 
are designed to alleviate the impact of spurious diffusion. 

In the second, "fluid-particle" method, the metallic- 
ity Z of a computational cell is calculated in the same 
way as in the first method, but the metal mass is cal- 
culated by summing the mass of the ejecta tracer par- 
ticles in the cell, and multiplying by the ejecta absolute 
metallicity ZgN =0.1. The tracer particle mass inside a 
cell is calculated with the cloud-in-cell method, with the 
weighting function taken to correspond to spreading the 
particle mass over a cubical kernel identical to the host 
cell, but centered on the particle. In the second method, 
only the cells with tracer particle presence contribute to 
the metallicity PDF. It should be kept in mind, however, 
that a metal particle displaced from the center of its host 
cell contributes mass to the neighboring cell that may be 
separated from the host cell by a shock transition or a 
contact discontinuity. 

In the third, "particle-particle" method, the metal- 
licity inside the cell is itself computed by summing up 

a step-function compositional discontinuity Z{x) = H{x) in an oth- 
erwise uniform medium traveling in the z-direction with velocity 
Vx > 0, and let Ax be the grid cell size. The code approximates the 
solution to the transport equation dZ/dt+vdZ/dx = by comput- 
ing cell boundary fluxes with the piecewise-parabolic method; here, 
we adopt a simpler scheme after assuming that just to the right of 
the discontinuity, the metallicity of every cell is much smaller than 
that of the cell immediately to the left. Let At denote the com- 
putational time step and S = vAt/Ax the dimensionless Courant 
parameter for fluid velocity, and assume that the discontinuity is 
slow, <5 <C 1. Then it is straightforward to show that after t/At 
steps, the metal-free region to the left of the traveling disconti- 
nuity will develop a spurious tail of nonzero metallicities of the 
asymptotic form 



\x/Axy 



(5) 



where the first factor is the binomial coefficient. The tail establishes 
itself with an unphysical speed ^ Ax/ At ^ v. 
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Figure 6. Metal (lower curve; Mz) and total baryon (upper 
curve; Mb) net mass flow rate through a sphere of radius 20 pc 
centered on the gravitational potential minimum. Dashed lines in- 
dicate outflows and solid lines inflows. The total mass of ejecta 
in the simulation is A/gN = 40 A/q and the metal mass is Mz = 
ZsnA'/sn = 4Af0. Net outflow reverses into an inflow earlier in 
the baryons because of the presence of cold fllaments delivering 
metal-free gas from the cosmic web into the halo center. 



the metal mass inside the cell (again, as in the second 
method, computed by adding up tracer particle mass in- 
side the cell and multiplying by .^sn), and dividing this 
by the sum of the metal mass and the fluid hydrogen 
and helium masses in the cell. Since the combined abun- 
dance of hydrogen and helium is near unity, the third 
method completely ignores the fluid metallicity as de- 
fined by the passive scalar advected on the computational 
grid. Ejecta particles as tracers of metallicity obviously 
do not suffer from the same spurious diffusion as the 
fluid metallicity, but they are affected by systematic is- 
sues of their own. In convergent flows such as near the 
insufficiently well resolved snowplow shell, metal parti- 
cles can get trapped in discontinuous fiow structures, 
which can produce a spurious increase of local metal- 
licity as computed by the third method. Thus, while 
the fluid metallicity is bounded from above by its max- 
imum value, .ZgN = 0.1, at the point of injection, the 
particle-based metallicity can exceed this value and even 
approach unity. 

In Figure [71 we show the met al-mass- weighted metal- 
licity PDF at four different times: first, in the Sedov- 
Taylor phase, and, then at approximately 1, 13, and 
40 Myr after the explosion. It is immediately clear that 
already at early times, the metallicity PDF exhibits a 
tail dMz/dhiZ ex Z extending toward very low metal- 
licities, where differences between the PDFs computed 
with the three methods became more severe. Spurious 
diffusion across compositional discontinuities in the fluid- 
based methods, and cloud-in-cell smearing in particle- 
based methods, both contribute to the tail. We also 
notice that the particle-particle method has developed 
a tail extending into the forbidden region, Z > Zg-^, a 
clear signature of spurious particle trapping. 



As the supernova remnant ages, metal mixing in the 
complex interior flow, which is numerical in the simula- 
tion but can be facilitated by turbulence in nature (e.g.. 
Pan et al. 2011), broadens the metallicity peak from its 
initial location at Z « Zgj^ to an approximate range 
10"^ < Z < Z„iax(i), where Z^^^{t) is a time- and PDF- 
extraction-method-dependent maximum absolute metal- 
licity. For fluid-based metallicity PDFs, we find that 
^max ~ (0.07, 0.02, 0.002) at t - igw « (1, 13, 40) Myr. 
We also notice a prominent narrow metallicity peak, built 
from cells belonging to the central, unresolved core that 
is accreting both from the pristine filaments of the cosmic 
web and from the ejecta-enriched supernova remnant fall- 
back. The peak is migrating slowly toward lower metal- 
licities, covering Zcorc ~ (2, 0.7, 0.3) x 10~^ in the same 
three epochs. The metallicity of the unresolved core is 
consistent with the average metallicity of the accreting 
fiuid given in equation (Hi). 

Finally, we turn to the energetics of the metal-enriched 
material. In Figure [8J we show the evolution of the total 
kinetic and internal (thermal) energy of fluid elements 
having absolute metallicity Z > 10~^. At Ikyr after the 
explosion, the kinetic energy decreases below its initial 
value when the ejecta start interacting with the ambient 
medium and the well-resolved reverse shock facilitates 
kinetic-to-internal conversion. Adiabatic and radiative 
cooling both contribute to the rapid decrease of inter- 
nal energy at 0.01 — 1.0 Myr after the explosion. The 
renewed rise of both energies in the old remnant reflect 
the fallback of ejecta-enriched fluid in the halo center. 

4. DISCUSSION 

4.1. Metal- Enriched Star Formation: Continuous, 
Bursty, or Self- Limiting? 

The prompt resumption of baryonic infall into the 
center of the cosmic minihalo suggests that star for- 
mation could resume on similarly short, ^ 1 — 5 Myr 
time scales. In the simulation presented here, the met- 
als that the supernova has synthesized begin returning 
to the halo center ~ 4 Myr after the explosion. The 
mean metallicity of the inflow, averaged over the pri- 
mordial streams and the metal- enriched collapsing su- 
pernova remnant streams converging in the center of the 
halo, is ^ 0.005 Zq assuming a net supernova metal 
yield of 4 Mq . If the metal-rich gas successfully mixes 
with the primordial gas, then, according to a prevailing 
belief, the metallicity of the combined streams is well 
sufficient to ensure that new star formation in this gas 
should be producing normal, low-metallicity Population 
II stars. The Population III to II transition cannot be 
ascertained in the present simulation because the simula- 
tion does not adequately resolve the cold gas at densities 
"■H ^ 100 cm^'^. This gas belongs to the central core and 
is colder than primordial gas at the same densities pre- 
ceding the insertion of the first star. It is at unresolved 
densities that the metallicity of the gas strongly affects 
the thermodynamic evolution a nd thus go verns th e out 
come o f fragme ntation (e.g. ,[Omukai et al. 2005; Santorol 



to] 

fc Shu ll 2006; Schneider et ar| PD06; Smith et al. 200^ 
Safranek-Shrader et al.||2010l assuming low mctallicities, 

z<omZQ). — 

The gaseous core in the center of the minihalo that is 
fed by metal-enriched streams can be compared to the 
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Figure 7. Absolute mctallicity distribution at four different times after the supernova. Tlie histograms show the fraction of total metal 
mass as a function of absolute metallicity calculated in three different ways: by computing the metallicity and metal mass from the fluid 
variables defined on the AMR grid (black curve), by computing metallicity from the fluid variables and metal mass from the Lagrangian 
particles tracing the supernova ejecta assuming ejecta metallicity of Zsn = 0.1 (red curve), and by computing both metallicity and metal 
mass with the Lagrangian particles (blue curve). The blue column in the first histogram shows the initial metallicity of the ejecta. The 
peak at Z ^ (2, 0.7, 0.3) X 10~* at t — isN = (1.187, 13.367, 40.127) Myr contains the fluid that has become incorporated in the central 
unresolved core (please see text for additional explanation). 



star- forming core of a giant molecular cloud (GMC). We 
attempt such comparison with an eye on the overall hy- 
drodynamics of the gas flow, putting aside any differences 
in chemical composition, dust content, CMB tempera- 
ture, etc., in the two systems. Dense GMC cores grow 
by accreting from larger-scale, lower-density cloud en- 
vironments which are supersonically turbulent. Smaller 
cores and their embedded forming star clusters merge 
with each other to form larger associations. Star for- 



mation in the cores is thought to be self-regulated, with 
the "feedback" from ongoing star formation controlling 
the amount of star-forming gas ava ilable and eventually 
leading to CMC dissolution (e.g., [Murray et al.||2010' 



a l][2010 

irpTTj ) 

H II re- 



Murray|[20TH [Goldbaum et al.|[20Tl[ [Lopez et al ||2Ulip 
I'he feedback entails, e.g., the pressure ot the H II re 
gions and the pressure of direct stellar radiation. While 
the baryonic mass of a cosmic minihalo is perhaps smaller 
than that of a GMC, the centrally concentrated under- 
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Figure 8. The kinetic (blue curve), internal (green curve), and 
combined kinetic and internal (red curve) energy of the metal en- 
riched fluid, defined as the fluid with absolute metallicity Z > 
lO"*^. 



lying dark matter halo renders the gravitational poten- 
tial depth in a minihalo's central gaseous core similar 
to that in a CMC's star-forming core. The underlying 
dark matter gravitational potential may make a cosmic 
minihalo less susceptible to star-formation-induced dis- 
solution than a CMC, assuming that the two systems are 
forming stars at the same rate. 

CMCs are thought to sustain star formation over tens 
of millions of years, in spite of the H II regions, su- 
pernovae, and other feedback activity that commences 
quickly following the formation of the first massive stars. 
The feedback's impact preceding the eventual dissolution 
of the cloud complex may be limited to reducing the av- 
erage star formation rate. We do not know whether star 
formation in metal-enriched minihalos and their descen- 
dent > lO"^ Mq halos can be sustained as in CMCs in 
spite of the feedback, or is the feedback from second- 
generation stars so effective as to evacuate baryons from 
the halo and prevent furt her star formatio n on cosmolog- 
ical time scales (see, e.g. 



Mori et alT2002|). It seems that 
d or 



these two scenarios can be evaluated only with attention 
to the detailed, small-scale mechanics of star formation 
and its feedback. A corollary of the first scenario would 
be the presence of multiple, chemically distinct stellar 
populations in the stellar system that forms this way. 

The statistics of low-luminosity stellar systems in the 
Local Croup in conjunction with theoretical modeling of 
structure formation in a ACDM universe suggests that 
star formation must have, to some extent, been sup- 
pressed in low-mass cosmic halos (e.g., Susa fc Umcmura 
2004al IRicotti & Cnedin"2005; Cnedin & Kravtsov 2006| 



Moore et al.||2006! 'Madau ct al. 20081 IBovill & Ricotti] 
2009, 2011; Koposov ct al. 2009; Muiioz et al. 2009; SaH 
vadori fc Ferrar a 2009; Busha et al.,2010, . Cao et al.,2010r 
Crifi:en et aTp OlO M accio et al.||2010| Kl^amoto et al 
20101 [Font etal.|2011i|Lunnan et al.|2012[ |Rashkov et al 



internally, as in CMCs, and externally, where the sup 
pression is driven by the heating of the ICM by the back 
ground ionizing radiation comi ng from oth er, more mas 
sive, less common halos (e.g., [Thoul fc W einberg 1996 
Kepner et al.||1997[ [Barkana fc Loeb 1999; Shapiro et] 



al.||20041 ISu sa fc Umemura 2004b; M esinger et al.. 2006. 
Mesinger & Dijkstra 2008; Okamoto et a,l.||2008p. The 
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heating could have shut off the inflow of fresh gas into 
the halo, eventually stunting star formation. The preva- 
lence of either suppression mechanism as a function of 
halo properties cannot be understood without first un- 
derstanding the detailed mechanics of star formation and 
feedback following the initial fallback of Population III 
supernova ejecta. 

4.2. The Implications of Uncertainties Related to the 
Nature of First Stars 

A source of uncertainty regarding the formation and 
nature of the first metal-enriched star clusters is the in- 
complete understanding of the precise character of metal- 
free Population III stars, and of the energetics and nu- 
cleosynthetic output of the associated supernovae. The- 
oretical investigations of Population III star formation 
have not yet converged and the precise functional form 
of the resulting initial mass function is not yet known. 
The evolution of metal-free stars is very sensitive to their 
poorly-understood internal rotational structure startin g 
with the protostellar phase (see, e.g., Stacy et al.||2011 ). 
Similarly, the explosion mechanism in metal-poor stars is 
not known. Each of the suggested mechanisms, including 
the standard delayed-neutrino mechanism hypothesized 
to drive Type II explosions, and the pair-instability and 
pulsational pair-instability mechanisms expected in rela- 



tively massive metal-free stars (s ee, e.g.,jQiatzopoulos & 
Wheeler||2012[ Yoon et aL||2012 and references therein), 
as well as the liiore exotic hypernova"-type mechanisms 
involving a (possibly magnetized) outflow from a central 
compact object, implies its own characteristic explosion 
energy and nucleosynthetic imprint. 

Many investigations of the formation of first stars 
and galaxies and their role in reionization have explored 
implications of the hypothesis that Population III stars 
had high masses, explosion energies, and metal yields. 
We find that the opposite limit of moderate masses, 
explosion energies, and metal yields might imply sharply 
divergent outcomes, not only for the parent cosmic mini- 
halos, but also for the overall pace of transformation of 
the early universe. A related issue is the nature of th e 
first galaxies (for a review, see |Bromm fc Yoshida|20TT ). 
Their assembly process is innuenced by the feedback 
from the Population III stars that formed in the minihalo 
progenitor systems, and this feedback in turn sensitively 
depends on the mass scale of the first stars. In the case 
of extremely massive stars, negative feedback is very 
disruptive, completely sterilizing the minihalos after 
only one episode of star formation. Only after ^ 10* yr 
did the hot gas sufficiently cool to enable the recollapse 
into a much more massive dark matter halo, where a 



2012 1. This suppression, however, could have arisen both 



second-gene ration starburst w ould be triggered (Wise & 
|Abel 2008b| [Creif et al.|[2010) ). This separation into two 
alstinct stages ot star formation would be smoothed out 
in the case studied here. Indeed, the first galaxies may 
then already have commenced their assembly in massive 
minihalos, as opposed to the atomic cooling halos (e.g.. 
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Oh fc Haiinan|2d02 ) implicated in the previous scenario. 



5. CONCLUSIONS 

We have carried out a cosmological hydrodynamical 
simulation designed to investigate the evolution of a cos- 
mic minihalo in the aftermath of the formation of the 
first, Population III star, assuming that the star has a 
moderate mass and explodes as a moderate-energy Type 
Il-like supernova, consistent with the recent downward 
revision of Population III star mas s estimates ( Stacy 
erari[20T0l [20121 IClark et al.|[20lT| |Greif et al.pJTTf 
Hosokawa et al.|2011[ ). We analyzed the dynamics of su- 
pernova ejecta and gas flow inside the minihalo. These 
can be compared to similar studies involving massive 
Population III stars exploding as ultra-energetic PISNe. 
Our main conclusions are as follows. 

The moderate-mass star and the moderate-energy su- 
pernova it produces inflict significantly smaller damage 
to the host minihalo than that inflicted by a massive 
star exploding as an ultra-energetic supernova. The star 
only partially photoionizes the host halo. The densest 
gas flowing from the filaments of the cosmic web into the 
star's parent cloud remains neutral and survives the su- 
pernova blastwave. This dense gas resumes accretion into 
the center of the minihalo only ~ 1 Myr after the explo- 
sion. After ~ 20 Myr, all spherically-averaged baryonic 
mass coordinates are moving inward and the accretion 
rate at the center of the halo is ^ 0.002 Mq yr~^. 

Following instability of the blastwave during the snow- 
plow phase, a fraction of the supernova ejecta starts 
falling back to the center of the halo. The ejecta fall- 
back accretion rate reaches a steady state after ^ 5 Myr. 
The accreting ejecta are incompletely mixed with the pri- 
mordial gas and are confined in thin sheets and filaments. 
Less than a half of the ejecta escape the virial radius; 
the escaping ejecta can be traced back to the outermost 
Rayleigh- Taylor fingers. 

The average metallicity of the accreting matter is in 
the range 0.001 — O.QIZq. This result depends on the 
assumed supernova energy and metal yield, but in this 
study, we did not pursue parameter space exploration. 
The metallicity of the accreted gas is sufficient to en- 
sure that new stars forming in the central core, unre- 
solved in the simulation, will have even lower character- 
istic masses. 

These results bring into focus metal-enriched star for- 
mation in cosmic minihalos and their immediate de- 
scendent halos. The character of post-supernova metal 
transport and gravitational fragmentation of the metal- 
enriched gas, and of the impact of subsequent star for- 
mation and supernovae, are both worth addressing in 
follow-up investigations. The early universe likely exhib- 
ited a range of different explosion settings, from mod- 
est to ultra-energetic. It will be intriguing to see how 
these high-redshift systems are matched with possible 
local fossil s, such a s the ult ra-faint d warf galaxie s (e.g. , 
Bovill & Ricotti|2009; Salvadori & Ferrara 2009, Tumlin 



Brodi 



le 



son 2010) and metal-poor globular clusters (e.g , 

S? Strader||2006| . From what we have already learned, 
Population 111 supernova feedback appears to play a key 
role in shaping early cosmic history and the nature of 
relic stellar systems. 
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